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iven  a hit  and  specific  target,  the  target's  vulnerability  depends  upon  the 
type  of  threat,  for  example,  kinetic  energy  (KE)  penetrator,  shaped-charge  jet 
or  high  explosive.  This  report  addresses  the  threats  associated  with  KE  rod 
penetrators  and  shaped-charge  jets.  The  questions  studied  are  penetration 
capability,  decreasing  penetrator  velocity  and  mass  through  the  armor,  hole- 
growth  rate  and  hole  size.  The  model  is  based  on  the  conservation  of  mass  and 
the  equations  of  motion  in  both  the  radial  and  axial  directions  and  is  used  to 
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predict  the  time  dependent  penetration  depth,  penetration  velocity,  radial 
velocity  imparted  to  the  target,  target  hole  profile,  and  the  erosion  of  the 
penetrating  projectile.  Incompressible  flow  is  assumed  for  the  axial  flow,  but 
target-strength  effects  and  viscous-shear  effects  are  included  in  the  analysis. 
Plugging  of  the  target  is  not  considered  in  the  present  model  end  all  target 
flow  is  assumed  to  occur  in  the  radial  direction.  This  radial  flow  (compress- 
ible but  uniform)  is  assumed  to  occur  in  three  different  regimes:  first,  target 
material  under  the  interface  flows  radially  outwards  (primary  regime),  then, 
eroded  penetrator  forces  additional  radial  flow  of  target  material  (secondary 
regime),  and  finally,  the  radial  inertia  induced  by  primary  and  secondary  flow 
is  absorbed  by  the  surrounding  target  material  (tertiary  regime) .-gjjue  to 
various  assumptions,  the  axial  and  radial  flows  are  uncoupled,  anothe  radial 
flow  can  be  obtained  once  the  penetration  and  penetrator  velocities  are 
determined.  Several  special  cases  of  axial  flow  are  considered  and  compare 
well  with  a variety  of  experimental  data.  Special  cases  of  radial  flow  will  be 
considered  in  a future  report. 
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INTRODUCTION 


Numerous  penetration  models  have  been  developed  in  the  past,  and  a 
large  number  of  these  and  the  corresponding  experimental  data  can  be 
foiond  in  References  1 through  14.  The  experimental  data  show  that 

^E.  M.  Pugh,  and  E.  L,  Fireman^  "Fundamentala  of  Jet  Penetration 
Carnegie  Institute  of  Technology^  Monthly  Report  OTB^4  (OSRD~4SS7) ^ 
Noveir^erj  1944. 

J.  Eiohelberger,  '^Re-examination  of  the  Theoriee  of  Jet  Formation 
and  Target  Penetration  by  Lined  Cavity  Chargee ^ " Carnegie  Institute  of 
Teohnologyt  CEL  Report  No.  2,  June^  19S4. 

^R.  DiPersiOi  J,  Simony  and  A.  MerendinOy  "Penetration  of  Shaped-Charge 
Je^  iT^^Mp^llia  Targets^"  BRL  Report  No.  1296,  September  1965. 

^F.  Allison  and  R.  Vitali,  "A  New  Method  of  Computing  Penetration 

-Shaped- C?zai*j70  Jets,"  BRL  Report  No.  1184,  January,  1963. 

^A.  Merendvno  and  R.  Vitali,  "The  Penetration  of  Shaped  Charge  Jets 
Into  Steel  and  Aluminum  Targets  of  Various  Strengths  (U),"  BRL  Memo- 
randum Report  No.  1932  (CLASSIFIED),  August,  1968.  (AD  #392672) 

^R.  DiPersio  and  J.  Simon,  "Resistanoe  of  Solid  Homogeneous  Targets  to 
Shaped  Charge  Jet  Penetration, " Ballistio  Research  Laboratories  Report 
No.  1417,  October  1968.  (AD  #841804) 

'^Huperveloaitu  Impact  Symposium  Proceedings,  Fourth  Symposium,  Vols.  I - 
III,  Air  Proving  Ground  Center,  Eglin  Avr  Force  Base,  September  1960. 

Defourneaux,  "Examination  of  the  Theoriee  Dealing  with  The  Penetra- 
tion of  A Projectile  into  a Plastic  Material,"  U.S.  Army  Foreign 
Science  and  Technology  Center,  FSTC-HT-23-1794-72,  13  March  1973. 

^D.  R.  Chidstman  and  J.  W.  Gehring,  "Semiannual  Report  on  Penetration 
MeaJumisms  of  High-Velocity  Projectiles,"  General  Motors  Report 
TR6 3-250,  31  December  1963. 

Johnson,  Impact  Strength  of  Materials,  Edward  Arnold,  London,  1972, 
Chapter  9. 

^^High-Velocity  Impact  Phenomena,  Edited  by  R.  Kinslow,  Academic  Press, 
1970,  Chap.  TV,  "Theory  of  Impact  on  Thin  Targets  and  Shields  and 
Correlation  With  Experiment,"  By  J.  W.  Gehring,  Jr.,  and  Chap.  IX 
"Engineering  Considerations  in  Hypervelooity  Impact,"  by  J.W.  Gehring, 
Jr. 

^^G.  Weihrauch,  "The  Behavior  of  Copper  Pins  Upon  Impacting  Various 
Materials  at  Velocities  of  50  - 1640  m/s,"  FSTC-HT-23-0004-74, 

10  September  1973. 

Tate,  "A  Theory  for  the  Deceleration  of  Long  Rods  After  Impact," 

J.  Mech.  Phys.  Solids,  Vol.  15,  1967. 

^‘*W.  A.  Allen,  E.  B.  Mayfield,  and  H.  L.  Morrison,  "Dynamics  of  a 
Projectile  Penetrating  Sand, " J.  of  Applied  Physics,  Vol.  28,  No.  3, 
March,  1957. 


strength  effects  are  important  for  both  KE  penetrators^*^  and  shaped- 
charge  jets^-»^»®  and  several  of  the  above  models  have  included  target- 
strength  effects.  Hovever,  the  majority  of  these  models  assume  that 
strength  effects  are  small  compared  to  the  effects  of  impact  velocity 
and  hence  consider  density  as  the  predominant  target  parameter.  Also, 
separate  models  have  been  developed  for  the  penetration  behavior  of 
shaped-charge  jets  versus  KE  penetrators  and  different  theories  exist 
for  the  two  types  of  penetrators.  Furthermore,  predictive  models  are 
lacking  for  the  mass  and  velocity  of  frontal  ejecta,  radial  hole  growth, 
hole  contour,  radial  velocity  induced  into  the  penetvator  at  the 
penetrator-target  interface  and  the  residual  mass  and  velocity  of  a rod 
penetrator.  The  model  introduced  in  this  report  is  applicable  to  both 
jets  and  rod  penetrators,  attempts  to  provide  predictive  models  which 
are  currently  lacking,  and  includes  the  effects  of  both  strength  and 
viscosity  of  the  target. 


II.  FORMULATION  AND  DERIVATION  OF  THE  MODEL 

The  intent  of  this  model  is  to  increase  the  accuracy  obtained  from 
earlier  models  (References  1 - 3,  7 and  9),  as  well  as  to  extend  the 
amouTit  of  information  derivable  from  simple  models.  In  an  earlier 
approach,  the  authors derived  an  approximate  model  which  relaxed  to 
known  solutions  for  shaped-charge  jets  and  penetrators,  and  included 
strength  and  viscosity  effects.  However,  the  radial  flow  of  target 
material  could  not  be  evaluated.  In  the  current  model,  both  axial  and 
radial  flow  of  the  target  are  considered  by  employing  conservation  of 
mass  and  momentum  sequentially  to  the  axial  and  radial  directions. 

There  are  several  ways  in  which  to  apply  the  conservation  laws  to 
the  penetration  problem.  In  the  differential  approach,  the  equations 
are  expressed  in  terms  of  unknown  stresses,  displacement  and  velocity 
fields.  This  would  yield  partial  differential  equations  analogous  to 
(although  possibly  simpler  than)  those  used  in  existing  hydrodynamic 
codes,  However,  the  usage  of  such  codes  for  a vulnerability  study 
would  be  prohibitive  in  both  time  and  cost.  Conversely,  a closed- form 
solution  would  be  readily  adaptable  to  a vulnerability  analysis,  and 
inspection  of  the  variables  contained  within  the  solution  would  provide 
insight  into  ways  of  either  enhancing  or  reducing  target  vulnerability 
(given  a hit).  Therefore,  in  an  attempt  to  obtain  a closed-form 
solution,  an  integral  approach  was  taken  to  the  problem. 

The  integral  approach  involves  approximating  the  real  problem  by  a 
series  of  regions,  in  each  of  which  the  internal  fields  are  assumed  to 
be  constant  and  whose  boundaries  are  acted  upon  by  some  average  surface- 
traction  at  the  centroid  of  the  surface.  These  regions  will  be  referred 

P,  Walters  and  J.  N.  Majerua,  "HyperoelooiUj  S.-rpaot  Models  for  Hole 
Growth  and  Geometry Prooeedin^s  of  the  Third  Annual  Vulnerability/ 
Survivability  Symposium  ADPAt  Naval  Amphibious  Base,  Coronada,  CA,  197?. 


8 


to  as  free-body  regions.  Applying  this  approach  to  the  penetrator 
problem,  the  mixing  of  the  penetrator  and  target  materials  in  the  inter- 
action region  and  the  convective  momentum-flux  terms  (spatial  gradients 
of  the  velocity)  are  ignored.  This  approach  serves  to  partially  decouple 
the  governing  equations  and  results  in  a system  of  ordinary  differential 
equations  with  respect  to  time. 

To  apply  this  approach  to  the  penetrator  problem,  we  first  need  to 
establish  a global  free-body  which  relates  the  penetrator/target  inter- 
face location  to  an  Eulerian  reference  plane.  This  establishes  the 
penetration  depth  P with  respect  to  the  top  of  the  target  and  defines 
the  axial  and  radial  velocities  as  dZ^^g/dt  and  dRgg/dt,  where  Rgg  and 
Zgg  are  the  center-of-mass  coordinates.  Figure  1 shows  the  penetrator 
at  three  global  stages:  initially  when  t » tjj  + At,  during  penetration 
(t  > tg  + At),  and  during  free-f light  at  the  same  time  t. 

The  top  of  Figure  1 shows  the  penetrator  at  some  time  At  after 
impacting  the  target  at  time  tQ.  During  this  time  At,  the  target 
r.aterial  under  the  penetrator  (to  depth  H)  is  replaced  by  penetrator 
material  which  has  slowed  down  to  a velocity  Uq.  This  region  of  pene- 
trator material  of  length  H and  of  penetrator/target  interface-velocity 
U is  referred  to  as  the  interaction  region.  The  details  of  the  forma- 
tion of  this  interaction  region  can  only  be  obtained  by  elaborate  hydro- 
dynamic-code calculations^®*^^  and  are  not  considered  in  the  present 
model . 

Note  that  the  time  At  is  assumed  to  be  small  compared  to  the  time  t 
associated  with  penetration  into  the  target,  i.e.,  t + At  “ t.  Also, 
the  interaction  height  H is  assumed  to  be  small  compared  to  either  the 


thickness  of  the  target  or  the  pen:jtrator- length  1^.  Flash  radiography  j 

and  hydrocode  calculations^ indicate  that  H is  on  the  order  of  the  5 

penetrator  radius  for  rod-type  penetrators-  Hence,  the  assumption  on  H 
should  be  valid  except  for  problems  involving  extremely  short  penetra-  ^ 

tors  or  thin  targets.  | 

The  second  and  third  global  stages  are  also  indicated  in  Figure  1.  j 

The  second  stage  refers  to  the  projectile  penetrating  the  material  ] 

whereas  the  third  stage  (free-flight)  is  merely  a reference  stage  ] 

required  for  stretching  shaped-charge  jets.  The  penetration  P is 

referenced  to  the  top  of  the  target  and  the  initial  penetrator  length  is  i 

(^o)-  ^ shaped  charge  jet,  Zq  (t^)  = Z^  (1  - Vj. (t^j/V^)  where 

Vq  is  the  tip  velocity,  Zq  is  the  distance  from  the  virtual  origin 


E.  Johnson^  "Three -DimenBiortat  Computc.tions  on  Penetrator^Target 
InteraationSt"  Ballistic  Research  Laboratory  Contract  Report  No.  338 ^ 
Api-il  1977.  (AD  #A041058) 

Kucher f "Computer  Study  of  the  Effects  of  Rod  Diameter  and  Target 
Thickness  on  the  Penetration  Process^  " Tech  Report  ARBRL-TR’-02046 ^ 
February,  1978. 
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PENETRATOR  AND  INTERACTION 
REGION  AT  TIME  {»o  ♦ At) 


PENETRATOR  AND  INTERACTION 
REGION  AT  TIME  t 


FREE  FLIGHT  PENETRATOR 
AT  TIME  t 


Figure  1.  Global  Penetrator  Stages 
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and  Vy(tp)  is  the  velocity  at  the  rear  of  the  jet  at  the  time  of  impact. 
Experiments  and  hydrocode  calculations  indicate  that  the  value  of  Vj. 
will  not  change  during  penetration  until  the  penetrator/targev  interface 
approaches  the  rear  of  the  jet.  For  a KE  penetrator,  Iq  ft.)  is  the 
initial  rod  length  and  V^ttg)  ■ at  the  time  of  impact.  Tf  the  rod 
is  decelerated  during  penetration,  the  vrlue  of  Vj.  will  change.  How- 
ever, Vp  will  be  approximately  equal  to  the  velc  Ity  V at  the  top  of 
the  Interaction  region. 

Based  upon  the  configuration  shown  in  Figure  1,  the  penetrator 
length  8,(t)  during  penetration  and  the  free-f light  penetrator  length 
tf(t)  at  the  same  time  are  respectively, 

“ ‘^o^^o  + At)  - AL  ♦ P .y*  V^(t)dt,  (1) 

t +At 
0 


^fCt)  = t^(t^  + At)  + (t  - (t^  + At))(V^  - V^(t^  + At),  (2) 


where  AL  is  the  length  of  penetrator  lost  during  the  formation  of  the 
interaction  region.  Also,  conservation  of  mass  of  the  free- flight 
penetrator  requires  that  the  rod  radius  Tj  be 


r?(t)  = ■ AL)/t£(t) 


(3) 


where  r is  rj  at  t = t^  + At  and  the  penetrator  is  assumed  to  have  the 
shape  ot  a circular  cylinder.  N< te  that  for  a KE  penetrator,  t£(t)  = 
S-oC^o)  - AL  and  r(t)  - r^Ct)  since  V^Ct^)  = for  free-f  light. 

Due  to  our  assumptions  concerning  At  and  H,  AL  should  be  small 
compared  to  Ji-gCto).  Hence,  tpCt^)  - AL  - At  ~ tg. 

The  penetrator' s velocity  V at  the  top  of  the  interaction  region  is 
taken  to  be  Vq  at  t » t^  + At,  and  the  penetration  depth  P is 
referenced  from  the  top  of  tho  interaction  region. 

Now  that  the  global  problem  has  been  formulated,  the  corresponding 
free-bodies  and  localized  flow  associated  with  the  interaction  region 
must  be  established.  The  flows  of  target  and  penetrator  materials  are 
envisioned  as  occurring  in  two  regimes.  The  first  regime  (regime  A) 
involves  a fixed  amount  of  penetrator  material  contained  within  the 
interaction  region,  moving  axially  with  velocity  U,  and  displacing 
target  material  associated  with  depth  H.  The  second  regime  (regime  B) 
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Involves  additional  penetrator  material  (due  to  the  relative  velocity 
V - U)  flowing  in  and  out  of  the  interaction  region.  The  penetrator 
flow  out  of  the  interaction  region  is  accomplished  by  displacing  target 
material  located  radially  around  the  penetrator.  Figure  2 shows 
schematic  diagrams  of  the  flow  associated  with  these  two  regimes  of 
flow. 

In  general,  part  of  the  target  material  under  the  interaction  region 
may  be  compressed,  and  the  remaining  material  flows  in  the  radial  and 
axial  directions.  The  compressed  target ^material  is  either  pushed  ahead 
of  the  penetrator  in  the  form  of  a so-called  shear  plug,  or  ejected 
from  the  back  of  the  target  in  the  form  of  debris.  The  target  material 
in  the  radial  region  around  the  penetrator  can  flow  in  both  the  axial 
and  radial  directions.  Similarly,  the  eroded  penetrator  mass 
(associated  with  the  flow  velocity  V - U)  can  flow  in  both  the  radial 
and  axial  directions.  However,  if  both  penetrator  and  target  materials 
flow  in  the  axial  and  radial  directions,  a large  number  of  coupled 
differential  equations  are  reqtured  for  solution.  Furthermore,  the 
partitioning  of  the  mass  flow  under  the  interaction  region  is  a complex 
phenomenon  and  may  involve  other  analyses.  Therefore,  for  the  present 
investigation,  the  following  assumptions  are  made  concerning  the  regime 
A and  B flows:  all  target  material  flows  radially  outwards  under  the 
interface  during  regime  A,  and  during  regime  B the  target  and  eroded 
penetrator  materials  flow  only  in  the  radial  direction. 

For  axial  flow  the  two  penetrator  regions  are  shown  in  Figure  3, 

In  the  present  approach,  the  values  of  the  penetrator' s initial 
velocities  (Vq  and  Vj.) , and  size  (r^^  and  , target  density  pene- 
trator density  Pg,  target  strength  Oj,  penetrator  strength  a,  and 
target  viscosity  pj,  and  interaction  viscosity  p are  assumed  to  be 
known.  The  unknowns  are  the  masses  M*  and  Mjj  of  the  two  regions,  the 
density  p of  the  penetrator  material  in  the  interaction  region,  the 
height  H of  the  interaction  region,  the  velocity  Vj  of  the  penetrator 
material  above  the  interaction  region,  the  velocity  V at  the  rear  of 
the  penetrator,  the  length  il  of  the  penetrator,  and  the  penetration 
depth  P.  The  penetration  velocity  U is  not  an  independent  variable 
since  it  is  related  to  penetration  through  the  definition  U = dP/dt. 
Also,  the  radius  r(t)  and  free-flight  length  )l£(t)  are  considered  known 
because  of  Equations  (2)  and  (3).  This  gives  eight  unknowns  and  one 
equation  (Equation  1).  Hence  seven  more  equations  are  needed. 

Because  of  our  definition  of  region  II,  ^ ® 

following  two  equations  apply: 

Mu  = ir  r^^Hp  = tt  r^(t^)  H(tQ)  (i) 


whereas  the  mass  of  region  I is 
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Figure  2.  Schematic  Diagrams  of  Penetrator  and  Target  Flow 
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Figure  3.  Free-Body  Configurations  for  Axial  Flow 


(5) 


Mj  - IT  r^(t)  pQ  A(t). 


Following  the  approach  of  Pugh,  Firemarf  and  Eichelberger^  for 
inertial  resistance,  equilibrium  of  axial  momentum^  for  region  II 
requires  that 


Mjj(U  + j H)  = IT  r^.^  (o  - o^)  + IT  r^.^p^  (Vj  - 
- IT  - IT  Hu^  2U 


(6) 


where  a superimposed  dot  denotes  differentiation  with  respect  to  time 
and  the  flow  shear- stresses  are  approximated  as  being  associated  with 
a linearized  Newtonian- fluid.  We  now  have  five  equations  for  the  eight 
unknowns.  Therefore,  three  additional  relationships  are  required. 


First,  the  height  H will  be  assumed  to  be  constant.  This  constant 
must  be  obtained  either  from  radiographic  observation,  or  hydrocode 
calculation^®* or  may  be  considered  in  a parametric  fashion.  Two 
additional  relationships  involve  two  sets  of  assumptions  made  separately 
for  jet  and  rod  penetrators,  i.e.. 

Vr  = V^,g(I)  (rod) 

(7a) 

V = V (jet) 
r ro 

(7b) 

and 

Vi  = Vj,q(I)  (rod) 

(8a) 

V,  = V + (V  - V (jet) 

I ro  o ro''  ' f 

(8b) 

where  VfjgCI^  = velocity  at  the  center  of  gravity  of  the  rod  penetrator 
and  VjQ  * Vj.(tQ).  The  justification  for  Equations  (7a),  (7b)  and  (8a) 
were  previously  given.  Equation  (8b)  corresponds  to  the  assumption  of 
a linear  velocity-gradient  in  the  jet.  Radiographic  observations^ 
show  that  this  is  a good  approximation  for  a variety  of  shaped-charges. 


4* 

If  not  all  of  the  eroded  penetrator  flowe  radially  during  regime  B, 
then  Mjj  changes  with  time  and  must  he  accounted  for  in  Equation  (6). 
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Equations  (7a)  and  (8a)  introduce  another  unknown  VQg(I)  which  can 
be  related  to  the  other  variables  through  the  conservation  of  axial 
momentum'*'  applied  to  region  I,  i.e.. 


. U)2  * o] 

- It  y P Vj,g(I). 


(9) 


Therefore,  the  axial  solution  for  shaped-charge  jets  involves  Equations 
(1)  through  (8b)  whereas  KE  penetrators  require  solving  Equations  (1) 
through  (9) . 


Note  that  the  height  of  the  interaction  region,  H,  need  not  always 
be  assumed  constant  but  could  vary  with  time.  One  applir ztion  would  be 
a penetrator  impacting  a rigid  (non-deforming)  target  where  H = H(t) 
but  the  penetration  is  zero  (see  Reference  18) . In  this  case  the 
problem  must  be  reformulated  slightly  with  region  II  being  attached  to 
the  top  of  the  target,  and  mass  Mjj  is  no  longer  constant. 

For  radial  flow  during  regimes  A and  B,  axial  views  of  the  regions 
of  interest  are  shown  in  Figures  4 and  5 respectively.  The  correspond- 
ing radial  target-regions  are  denoted  as  regions  III  and  IV  respectively. 
These  regions  are  of  height  H and  located  at  depth  P within  the  target. 
The  penetrator/target  interface  reaches  this  depth  at  some  time  tp. 

The  radial  extent  of  regions  III  and  IV  are  denoted  as  Rp*  and 
Rp’^*  respectively. 

During  flow  regime  A,  (see  Figure  4),  target  material  flows  out 
from  under  the  interface  and  compresses  the  target  material  contained 
within  region  III.  The  density  p^*  changes  with  time  but  is  assumed  to 
be  uniform  throughout  region  III.  During  flow  regime  A,  the  radial 
boundary  R_,*  is  assumed  to  be  stationary.  This  flow  regime  lasts  for 
time  t*  = H/U. 


During  flow  regime  B (see  Figure  5),  the  eroded  penetrator  flows 
radially  outwards  and  compresses  the  target  material  in  region  IV. 

The  corresponding  density  Pj**  changes  with  time  but  is  assumed  to  be 
iiniform.  The  radial  boundary  Rp**  is  assumed  to  be  stationary  and  flow 
regime  B lasts  for  a time  t*. 


'^Appendia  I lists  a second  approach  for  region  I which  is  based  upon  a 
global  momentum  conservation  equation  where  region  II  is  viewed  as  a 
shock  front  with  H 0. 

Vitkins  and  M.  Guinam,  "Impact  of  Cylinders  on  a Rigid  Boundary t" 
J,  of  Appl,  Fhys.t  ^ol.  44,  No.  3,  March  1973. 


Figure  4.  Two  Axial  Views  of  Radial  Seginent  for  Flow  Regime  A 


There  are  six  unknowns  associated  with  each  regime  of  flow,  i.e., 
the  regions  III  and  IV.  These  unknowns  are  the  following;  the  radial 
locations  r*  and  r**  of  the  penetrator/target  interface,  the  radial 
locations  Rq*  and  Rq**  of  the  boundaries,  the  locations  R*^  and  R*g  and 

velocities  V*^  and  V**  of  the  cent er-of -mass  of  the  regions,  and  the 

masses  M*  and  M**  and  densities  p*  and  pij*  of  the  target  material. 

Therefore,  a total  of  12  equations  are  required  to  solve  for  the  regime 
A and  B flows. 

Eight  equations  can  be  obtained  from  the  geometry  and  definition  of 
regions  III  and  IV, 


segment 

“"lunlt  = " ClOb) 

segment 


V*  = R* 
CG  ‘^CG 


V**  = R** 
^CG  CG 


R*  = — 
CG  3 


2 + r*  R^*  + (r*)^] 


(R  * + r*) 
o 


+ T* 


(R** 
^ o 


+ r**) 


(lla) 

(lib) 

(12a) 

(12b) 

(13a) 

(13b) 


The  equations  for  r*  and  r**  stem  from  equating  the  axial  penetrator 
flux  into  the  reference  regions  to  the  psuedo  radial -fluxes  (see  Figures 
4 and  5).  This  gives  the  following  equations,  for  regions  III  and  IV 
respectively. 
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(rod) 


(14a) 


r*  yet  - t ) U r^^/H 


Vr  ^ H U (V  - V )t  1 

o 0 o o ro-*  f... 

H(V  - V ) I + (V  - V )t 
0 vo'^  L o '■  o ro"^  p. 


r**  = y (t  - t ) (V  - U)  r^^/H  (rod) 


(14b) 


(15a) 


4. 


i (V  - U) 


An 


ri  + (V  - V )t  1 
0 o ro"^ 

A + (V  -V  )t' 

_ o o ro"^  pl_ 


(jet) 


(15b) 


where  tpj  = tp  + t*,  Equations  (14a)  and  (14b)  are  limited  to  the  time 

interval  t^<  t < + t*,  and  Equations  (15a)  and  (15b)  are  limited  to 

the  interval  t , < t < t , + t*. 

pi  pi 

The  last  two  equations  deal  with  the  conservation  of  radial 
momentum  at  various  times  during  the  flow  regimes,  For  flow  regime  A, 
the  instantaneous  change  of  axial  flow  to  radial  flow  at  time  t = tp 
introduces  a mathematical  singularity  in  the  radial  momentum  equation, 
i.e. , lim  f*  •♦  ».  Therefore,  for  flow  regime  A,  conservation  of 

t-»-tp 

radial  momentxim  will  be  applied  only  at  the  end  of  this  flow  regime, 
i.e.,  at  t = t_  + t*.  Again,  following  the  Pugh-Fireman-Eichelberger 
approach  for  inertial  resistance,  the  balance  of  radial  momentum  yields, 

M*  ^Cg|  •{■'';«  * (''cg’'  '>To1 

t = t + t*  ' 

P 


- [(r*)2  - (r*)2]  yyn  (16) 


+ r*  H [o  + (V*^  - f*g)2pj| 

t = tp  + t* 
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where  > 2/3  f*  and  the  material's  resistance  is  represented  by  a 
strength  acting  normal  to  the  surfaces,  and  a linearized  Newtonian 
shear-stress  acting  tangentially  to  the  surfaces.  Note  that,  because 
of  Equations  (10)  through  (IS) , the  only  unknown  in  Equation  (16)  is 
the  radial  location  RJ  which  is  constant  at  time  t « tp  + t*.  There- 
fore, Equation  (16)  represents  an  algebraic  equation  (hon- linear)  for 
obtaining  R*. 

The  conser  cion  of  radial  raomentiam  is  applied  to  region  IV  so  that 
the  change  in  momentum  between  times  tp  + t*  and  tp  + 2t*  equals  the 
external  forces  acting  over  this  time  interval,  i.e., 

+ 2t* 

d(MV)  = M**  V**  (tp  + 2t*)  - M*  (tp  + t*) 


t + t* 

P 

- [(R**)^  - (r**)^]  v**/H 

+ r**  H [0  + (V**  - h**)^  p^]|  dt 


where  ffg  = 2/3  f**.  Again,  Equations  (10)  through  (15)  yield  the  time 
dependency  for  the  integrand  of  Equation  (IT),  and  the  solution  to 
Equation  (16)  yields  the  value  of  M*  and  VCG-  Hence,  Equation  (17)  is 
an  algebraic  equation  for  the  unknown  constant  R**. 

Figure  6 shows  radial  segments  of  target  and  penetrator  materials 
as  viewed  at  three  different  times.  The  top  sketch  refers  to  the  target 
and  penetrator  materials  at  the  end  of  flow  regime  A (t  = t..  + t*) ; 
the  middle  sketch  corresponds  to  the  materials  at  the  end  or  flow 
regime  B (t  = tp  « 2t*) ; and  the  bottom  sketch  views  the  materials  at 
some  time  after  the  interaction  region  has  passed.  This  latter  stage 
of  radial  flow  will  be  referred  to  as  flow  regime  C.  During  this  flow 
regime  the  outward  flow  for  a unit  segment  of  fixed  target  and  pene- 
trator masses  are  impeded  by  the  target  material  remaining  in  the 
radial  direction.  The  following  11  unknowns  are  associated  with  flow 
regime  C:  the  constant  masses  Mp  and  NU.  and  time-dependent  densities 
p and  p„  of  the  penetrator  and  target  materials,  the  radial  location 
rp  of  the  penetrator  material,  the  radial  location  Rj  of  the  penetrator/ 
target  interface,  the  radial  location  R of  the  interface  between  the 
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moving  and  stationaiy  target  material,  and  lastly,  the  radial  locations 
Rqq  and  Tqq,  and  respective  velocities  Vj  and  Vp,  of  the  center^of- 

masses. 

Based  \ipon  the  geometry  and  definition  of  flow  regime  C.  the  follow- 
ing 8 equations  can  be  applied. 


Hf  ■ h Rto/2 

(18a) 

Mp  - (r**  (tp  ♦ 2t*)fH  p^/2 

(18b) 

H [R^  - Rj^]  - Wj. 

(19a) 

H 0 [R,^  - - Mp 

(19b) 

2 [R^  * RRj  * Ri'l 
‘^CG  3 [ R + Rj  J 

(20a) 

- + r R-  + r ^"1 

r = 2 J 2J 2_ 

^CG  3 L Rj^rp  J 

(20b) 

* '^CG 

(21a) 

'^p  “ ^CG  • 

(21b) 

Equilibrium  of  radial  momentum  for  the  moving  target  and  penetrator 
segments  yields  the  following  two  eqiiations. 

“r  \ " H * ‘’to  V’ 

- Uj  (R^  - Rj2)/H  (22) 

+ Rj  H [a  + p (V^  - Vp)2] 
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and 


M V 
P P 


- RjH 


[Oy  + Pj 


(V, 


T ■ V ^ 


(23) 


«We  have  10  equations  and  11  unknowns.  Therefore  an  additional 
assumption  is  needed.  Depending  upon  the  penetrator,  the  authors 
propose  either  of  the  following  two;  (a)  assume  rp  = 0 and  hence  the 
penetrator  material  is  expanding  into  the  hole  or;  (b)  assume  P = Pq 
and  hence  the  radial  shape  of  the  penetrator  changes  as  the  hole  expands. 
Based, upon  either  assumption,  the  location  and  velocity  of  the  penetra- 
tor segment  is  directly  related  to  the  interface  location  Kj.  Assump- 
tion (a)  relates  to  Rj  via  Equations  (20b)  and  (21b),  whereas  assumption 
(b)  relates  to  Rj  through  Equations  (18b),  (19b) , and  (21b) . 


Equations  (18)  through  (21)  and  either  of  the  above  assumptions  can 
be  combined  with  Equations  (22)  and  (23)  to  yield  two  coupled  differ- 
ential equations  for  the  unknowns  Rj  and  R.  The  initial  conditions 
associated  with  these  differential  equations  come  from  the  solution  of 
Equations  (10)  through  (17),  i.e.. 


R (tp  + 2t*)  * R** 

R,  (t  + 2t*)  r**  (t  + 2t*) 

F r 

Vt  (S  " “ ''CG  ^"p  " 

Vp  (tp  . 2t*)  = f-  (tp  . 2t*) . 


(24) 


III.  DISCUSSION 


Due  to  the  assumed  conversion  between  axial  and  radial  flows,  the 
axial-flow  equations  are  independent  of  the  radial-flow  equations. 
Therefore,  the  hierarchy  of  the  solvition  scheme  for  the  above  equations 
is  as  follows.  First,  for  the  jet,  one  must  solve  for  the  penetration 
velocity  U from  Equation  (6)  (with  dH/dt  = d^H/dt^  = 0)  with  Equations 
(1)  and  (8)  being  used  for  the  jet  length  I and  penetrator  velocity, 
whereas  for  a rod.  Equations  (6)  and  (9)  must  be  solved  simultaneously 
for  U and  V,  with  Equations  (1)  and  (5)  being  used  for  I and  Mj.  Once 
U and  V are  determined.  Equations  (14)  and  (15)  can  be  evaluated  in 
terms  of  a penetration  time  t-,  and  algebraic  Equations  (16)  and  (17) 
yield  the  unknown  radial  boundaries  Rj  and  Rj*.  With  Rq  and  Rq*  known, 
Equations  (18)  through  (24)  can  be  evaluated.  Then,  Equations  (22)  and 
(23)  must  be  solved  simultaneously  for  R and  Rj.  Note  that  Rj 
corresponds  to  the  instantaneous  hole  profile  located  at  penetration 
depth  P.  The  solutions  for  R and  Rj  depend  implicitly  upon  the  time  tp 
and  usage  of  the  known  U (=  dP/dt)  permits  the  corresponding  depth  of  ^ 
penetration  to  be  calculated.  Hence,  one  obtains  the  hole  growth  at 
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depths  of  penetration  associated  with  preselected  values  of  tp. 

Because  of  the  non- linearities  of  the  above  set  of  equations,  a 
general  closed-form  solution  for  either  a rod  or  jet  was  not  evident  to 
the  authors.  For  axial  flow  an  elaborate  Integral-differential  equation 
for  V Is  obtained  for  a XE  rod,  and  a second-order  differential  equation 
for  P develops  for  a shaped-charge  jet.  The  second  order  equation  Is 
extremely  cumbersome  since  the  radius  r Is  a function  of  time  and  the 
time  derivatives  of  momentum  involve  fj(t),  i.e.. 


i,  (t) 


- r V 
0 o 


* * Vh'i 


372 


’ * t/t 


C25) 


where  the  definitions  of  and  t^  have  been  substituted  into  Equation 
(25).  Interestingly,  Equation  (25)  shows  that  the  time  variation  {"(t) 
will  be  small  for  the  condition  of  a large  value  of  t^,  i.e.,  large 
standoff-distances  2^  or  slow  velocities  V^.  For  this  condition,  the 
jet’s  radius  can  be  considered  to  be  approximately  constant.  This 
greatly  simplifies  the  conservation  equations  and  enables  approximate 
closed-form  solutions  to  be  obtained.  In  fact,  numerical  solutions  of 
Equations  (6),  and  (9)  are  underway  using  the  method  of  weighted 
residuals^®. 

Whereas  the  equations  for  radial  target  flow  (Equations  (23)  and 
(24))  are  not  complicated  by  the  expressions  for  )!,(t)  and  r(t),  they 
are  second-order  with  respect  to  both  R and  Rj  and  are  extremely  non- 
linear because  of  the  relationships  of  Equations  (20)  and  (21).  Again, 
no  closed-form  solution  was  obvious  to  the  authors.  An  order-of- 
magnitude  analysis  and  some  special-case  solutions  are  currently  being 
investigated  and  will  be  reported  at  a later  date. 

For  the  case  of  axial  flow,  special  closed-form  solutions  can  be 
obtained  if  other  assumptions  are  made.  Appendix  I gives  the  axial - 
momentum  equation  for  a KE  penetrator  when  region  II  is  viewed  as  a 
shock  interface.  Applying  the  jump  conditions^®  to  this  interface, 
one  obtains  a relationship  between  V and  U.  Then,  assuming  a constant- 
length  penetrator,  the  penetratioi.  and  center-of-mass  velocities  are 
obtained  as  a function  of  time  for  inviscid  flow.  The  influence  of 
viscous  effects  is  shown  via  an  approximate  solution. 


A.  Finlayaon  and  L,  E.  Seriven^  "The  Method  of  Weighted  Reeiduala  ~ 
A BevieWj"  Applied  Meahaniaa  Revieway  Vol.  19y  Wo.  9,  September*  1966. 

CviateaaUy  "Dynamia  Plaatioityy " Chapter*  Vllly  Vol.  4,  Wor*th- 
Holland  Publiahing  Companyy  Am8ter*dmy  196?. 
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Appendix  II  gives  the  solution  for  a steady- state  case  and  one 
obtains  the  familiar  shaped-charge  jet  relationships.  Also,  an  order- 
of -magnitude  analysis,  based  on  the  assumption  that  the  normalized  jet 
radius  is  small,  reveals  that  the  familiar  shaped-charge  jet  relation- 
ships can  be  obtained  without  assuming  axial  steady-state.  The  order- 
of-magnitude  analysis  also  shows  the  relaxation  to  the  Allison-Vitali** 
and  DSM^  shaped- charge  jet  theories.  Appendix  III  contains  the  solutions 
to  the  special  cases  of  Johnson's  KE  penetrator^®  and  constant  velocity 
penetrators  which  merely  erode  during  penetration.  Appendix  IV  gives 
the  general  solution  (with  strength  and  viscosity)  to  the  problem  of  a 
non-eroding  penetrator  which  is  decelerated  by  the  target,  e.g.,  a 
steel  bullet  into  sand  or  clay. 

It  should  be  noted  that,  other  than  the  steady-statr  solution,  all 
the  special  solutions  exhibit  penetration  versus  time  behavior  analogous 
to  that  observed  for  a shaped-charge  jet  (see  Figure  7) . The  data  in 
Figure  7 have  been  normalized  using  Pmax  where  t/tg  = 1 corre- 

sponds to  the  approximate  time  when  tne  jet  starts  to  break  up  into 
discrete  segments.  Note  that  the  penetration  proceeds  in  a continuous 
fashion  throughout  the  entire  time  domain.  Also,  for  moderate  changes 
in  standoff  condition,  there  appears  to  be  only  a slight  change  in  the 
penetration  verstis  time  behavior  (see  Figure  8).  Therefore,  even  though 
the  model  was  developed  for  a continuous  jet,  it  is  the  authors  con- 
tention that  the  model  should  be  applicable  to  a particulated  jet  as 
long  as  the  spacing  between  jet  particles  is  small  compared  to  the 
particle  lengths.  When  this  is  not  the  case,  one  can  consecutively 
apply  the  KE  rod  equations  to  each  particle  in  the  jet. 

In  general,  all  of  the  solutions  yield  time-dependent  functions. 
However,  maximum  penetration  Pj„  can  be  obtained  by  evaluating 
P(c  ■+•  t^)  where  tj^  is  the  time  where  the  penetration  process  stops. 

The  penetration  procesj  stops  when  the  penetration  velocity  U goes  to 
zero.  Therefore,  by  equating  the  expression  for  U equal  to  zero,  one 
can  solve  for  t^(l).  Another  value  of  maximum  time  tn,(2)  can  be 
obtained  by  equating  the  penetrator  length  (Equation  (1))  equal  to  zero. 
Note  that,  depending  upon  target  and  penetrator  conditions,  both  values 
of  tj^  might  approach  infinity.  By  comparing  these  values  t^,,  insight 
into  the  penetration  process  can  be  obtained.  If  t,„(l)  - tjn(2) , this 
implies  that  the  penetrator  is  completely  used  up  during  penetration 
and  the  penetration  process  stops  when  the  penetrator  is  dissipated, 
i.e.,  no  secondary  penetration  occurs.  If  t,n(l)  > t^(2),  penetration 
continues  after  the  penetrator  is  dissipated  and  secondary  penetration 
occurs,  e.g.,  hypervelocity  impact  into  a lead  target.  Lastly,  if 
t,n(l)  < ^m(2)  the  penetration  process  stops  before  the  penetrator  is 
completely  expended,  e.g.,  penetration  by  a shaped-charge  jet,  or  low 
velocity  impact  by  a tungsten  rod  into  an  aluminum  target. 

Finally,  in  Appendices  II  through  IV  the  various  special-case 
closed-form  solutions  are  compared  to  the  experimental  data  from 
References  10  and  12  through  14.  In  all  cases  for  rods,  the  agreement 
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Figure  7.  Experimental  Penetration  versus  Time  for  a Copper 
Shaped-Charge  Jet  into  Various  Target  Materials 


between  the  special  case  solutions  and  the  experimental  data  is  good. 
Note  that  the  values  for  strength  and  viscosity  were  obtained  a priori 
from  tabulated  values.  No  attempt  was  made  to  curve-fit  the  data  by 
changing  either  c or  y. 

For  the  case  of  the  shaped-charge  jets,  the  authors  were  interested 
in  showing  the  effects  of  various  strengths  and  viscosities  upon  the 
predicted  results  for  the  approximate  models.  The  approximate  models 
were  also  judged  according  to  how  closely  they  matched  the  experimental 
data.  Both  the  steady-state  and  the  constant -velocity  models  predict 
penetrations  considerably  less  than  measured  when  the  rod  length  is  set 
equal  to  the  initial  standoff  2q.  If  the  fully  stretched  jet-length 
(■^bC'^o  ■ '^ro^^  used,  the  predicted  penetrations  are  approximately 
equal  to  the  observed  penetrations,  except  for  aluminum.  However,  using 
this  length  gives  the  erroneous^* prediction  that  the  jet's  penetra- 
tion is  independent  of  the  standoff  distance  Zq.  The  constant -length 
model  gives  results  in  fair  agreement  with  the  experimental  results 
when  Iq  = Zq.  This  model  indicated  that,  of  the  four  target  materials, 
the  penetration  calculations  for  only  aluminum  and  lead  were  sensitive 
to  the  target  strength  (over  the  range  selected).  Conversely,  the 
predicted  results  for  steel  are  highly  dependent  upon  the  viscosity 
value  whereas  aluminum  is  only  slightly  sensitive  to  the  chosen  value. 
Table  IV-1  lists  the  penetration  values  discussed  above. 


IV.  SUMMARY  AND  CONCLUSIONS 

A penetration-hole  growth  model  was  developed  for  both  KE  rods  and 
shaped-charge  jet  penetrators.  This  model  involves  the  penetrator,  an 
interaction  region  located  above  the  target/penetrator  interface,  and 
three  radial-flow  regimes,  and  includes  the  effects  of  stretching  jets, 
material  strengths  and  viscous  shear.  The  interaction  region  is  assumed 
to  form  during  some  small  time  t*  after  impact  and  the  axial  flow  of 
the  target  and  penetrator  materials  is  assumed  to  be  incompressible. 

The  radial  flow  of  target  material  is  uncoupled  from  the  axial  flow  of 
penetrator,  i.e.,  the  solution  for  axial  flow  (penetration)  is  indepen- 
dent of  the  radial  flow.  Radfal  flow  is  assumed  to  take  place  in  three 
regimes.  First,  the  target  ma‘,erial  under  the  interface  flows  radially 
outwards  and  is  compressed  (pri7i..iry  flow  regime) . Next , the  eroded 
penetrator  material  (due  to  the  velocity  difference  Vj  - U)  causes 
additional  radial  flow  and  compression  of  the  target  material  (secondary 
flow  regime).  Lastly,  the  radial  target  momentum  induced  during  the 
primary  and  secondary  flow  regimes  is  absorbed  by  the  surrounding  target 
material  (tertiary  flow  regime) . 

The  normalized  form  of  the  governing  equations  for  axial  flow  indi- 
cate that  the  appropriate  non-dimensional  groups  are  the  density  ratio, 
the  Best  number,  and  the  Reynolds  number.  Due  to  the  non-linearity  of 
the  equations,  the  general  case  must  be  treated  numerically.  However, 
several  cases  are  presented  which  yield  exact  solutions. 
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The  first  special  case  is  axial  steady-state  penetration  with 
negligible  strength  and  viscous  effects.  This  case  corresponds  to  the 
classical  shaped-charge  jet  approach  and  the  model  yields  both 
Bernoulli's  equation  and  a maximum  penetration  which  are  in  agreement 
with  the  classical  results^"'*. 

Also,  an  order- of-magnitude  analysis  reveals  that  Bernoulli's 
equation  represents  the  dominant  order  terms  in  the  axial  momentum 
equation  when  the  ordering  analysis  assumes  that  the  normalized  pene- 
trator-radius  is  small.  Also,  the  Allison-Vitali**  continuous  jet  theory 
follows  from  the  order-of-magnitude  analysis. 

The  second  special  case  considers  Johnson's  rod  penetrator^*^  (a 
penetrator  of  constant  mass).  The  model's  results  are  analogous  to 
that  presented  in  Reference  10  for  inviscid  flow  and  negligible 
strength  effects. 

The  third  special  case  deals  with  a uniform  velocity  penetrator 
which  does  not  decelerate  but  does  erode  away  during  the  penetration 
process.  The  predicted  penetration  appears  to  follow  the  general  trend 
of  experimental  data,  i.e.,  the  slope  of  the  penetration  versus  time 
behavior  is  monotonically  decreasi^n?;.  In  fact.  Figures  9 and  10  indi- 
cate good  agreement  between  the  theory  and  the  experimental  results 
from  Reference  13. 

It  should  be  noted  that  all  of  the  above  special  cases  could  be 
readily  extended  to  allow  the  penetrator  radius  to  vary  witli  time,  and 
to  include  both  strength  and  viscous  effects. 

The  last  special  case  considers  a rigid  rod  penetrator  which  is 
decelerated  by  the  target  but  does  not  erode  away.  The  model  shows  that 
the  penetrator  velocity  V must  equal  the  penetration  velocity  U for  this 
special  case.  In  this  case,  a closed-form  solution  is  obtained  for 
U(t),  penetration  P(t),  time  for  penetr'Jtion  t_(l)  , and  maximum  penetra- 
tion P,„.  All  of  these  solutions  include  the  effects  of  both  shear 
viscosity  and  material  strength.  The  penetration  vs  time  behavior 
follows  the  experimental  trends  and  shows  that  the  slope  decreases  with 
increasing  target  strength  and  viscosity. 

The  results  of  this  case  (Appendix  * ..s  shown  in  Figures  11 

through  13,  where  again  the  agreement  with  the  experimental  data  from 
References  14  and  12  is  good  and  Figure  12  matches  the  trend  of  the 
shaped-charge  jet  data  given  in  Figures  7 and  8. 

In  summary,  it  appears  that  the  model  is  capable  of  yielding  known 
results  for  several  special  cases,  These'^sults  imply  the  generality  of 
the  model.  In  addition  to  predicting  Pjj,  the  model  predicts  (requires 
an  approximate  method  of  solution  for  the  general  case)  the  penetration, 
penetrator  erosion,  penetrator  mass  and  velocity,  penetration  velocity, 
and  radial  hole  size  and  growth  rate,  all  as  a function  of  time. 
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Figure  10.  Comparison  of  the  Experimental  Penetration  Data 
of  [13]  with  Appendix  III  Results 
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from  [14]  with  Appendix  IV 
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Figure  13.  Comparison  of  the  Experimental  Penetration  Data 
from  [12]  with  Appendix  IV  Results 
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APPENDIX  I 


AXIAL-MOMENTUM  EQUATION  FOR  A ROD  WHEN  REGION  II 
IS  VIEWED  AS  A SHOCK  INTERFACE 

In  this  approach,  the  axial -momentum  equation  is  of  the  following 
foim, 


^(Mj  V(.gCi) 


(i-n 


where  the  difference  in  surface  tractions  (oj , 02)  across  the  shock 
interface  (moving  with  velocity  U with  respect  to  the  Eulerian  reference) 
must  be  obtained  from  the  jump  conditions'^.  The  flow  shear-stresses 
acting  along  the  embedded  length  of  penetrator  is  represented  by  a 
linearized  Newtonian- fluid.  The  jump  conditions  associated  with 
continuity  and  momentum  yield  the  following  two  equations^^. 


Cl  - 02 


Pi  Pj  Vj  + P2  P2  V2 


(1-2) 


?!  Pj  = P2  P2  (1-5) 

where  Pj  is  velocity  of  propagation  of  the  interface  with  respect  to 
particles  in  region  1 (penetrator),  i.e.,  - U + Vqq(I),  P2  is  the 

velocity  of  propagation  of  the  interface  with  respect  to  particles  in 
region  2 (target),  i.e.,  P2  * U;  is  the  particle  velocity  in  region 

1,  i.e.,  Vj  » V^q(I);  and  ^2  particle  velocity  in  region  2,  i.e., 

V2  = 0.  Therefore,  Equations  (1-2)  and  (1-3)  become, 

^ - “2  ' “l  ''cg”’  - ''cg''” 

Pj  CV(.5C1)  - U)  . Pj  U.  (I-S) 

Equation  (1-5)  can  be  solved  for  V^q(I)  to  give, 

Vj,q(I)  = (1  + P2/P1)  U (1-6) 

and  this  relationship  and  Equation  (1-5)  substituted  into  Equation 
(I-l)  yields. 
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(1-7) 


^(Mj  V^(l))  . - -rj*  Pj  vJgCD/Cl  * pj/B,) 

- ir  y P Vj,j,(I) . 

Note  that,  if  the  densities  of  the  target  and  penetrator  materials 
or  either  side  of  the  shock  are  taken  to  be  the  initial  material 
densities,  then  V^q(I)  ■ [1  + V^l  U for  the  shock  relationship  whereas 

Vj  » (1  + Y)  U for  the  steady-state  jet  expression  (see  Appendix  II). 

e 

Now,  if  the  rod  length  is  assumed  to  be  constant  (Mj  « 0),  several 

approximate  solutions  can  be  obtained.  First,  if  the  viscous  effects 
in  the  interaction  region  are  negligible,  Equation  (1-7)  can  be  solved 
for  Vj.q(I),  i.e. , 


V i 
o o 


V^(t  - t^)/(l  + Y^) 


(1-8) 


Equation  (1-8)  is  now  substituted  into  Equation  (1-6)  and  integrated 
to  give  the  following  penetration  versus  time  relationship. 


P = tn  [1  + Y . 

i ■ 


(1-9) 


Note  that  there  is  no  finite  value  of  P^g^.  associated  with  this 
inviscid  model. 

When  the  viscosity  is  not  negligible,  one  can  transform  (letting 
P = X(P))  Equation  (1-7)  and  obtain  the  following  first  order  equation, 


Ae"^^  /^o  +B  - CP  = D P 


(I-IO) 


where 


* = \ “o  - ‘o 

B = t V 
o 

C * v 

D B Y^  r ^ 
o 


and  is  the  initial  interface-velocity  (the  particle  velocity 
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associated  with  the  striking  velocity  and  the  Itugonlot  behavior  of  the 
target  and  penetrator  materials).  An  explicit  solution  of  Equation 
(I-IO)  was  not  evident  to  the  authors.  However,  if  the  exponential  is 
expanded  and  the  first  three  terms  retained  (good  approximation 
P Y^/*-o  ^ following  differentipl  equation  is  obtained, 


A + B - P 


D P 


(I-ll) 


Equation  (I-ll)  can  be  viewed  as  a quadratic  for  P^,^  based  on  the 
penetration  velocity  (P)  going  to  zero.  The  corresponding  solution  is, 

[^^4=35]  (.-.2, 

2 A 

where  A 5 (1  - v Hq/Tq  Uq)  and  Equation  (1-12)  is  subject  to  the 
restriction  that  PY^/*'q  < 1-  Equation  (1-12)  yields  non-imaginary 
solutions  and  satisfies  the  restriction  when  A < 0,  i.e., 

V a > Yi^  r ^ U . (1-13) 

0 0 0 


Furthermore,  since  P^gy  must  be  a positive  number,  the  negative  sign 
in  Equation  (1-12)  must  be  selected. 

Equation  (I-ll)  can  also  be  rearranged  to  give  the  follovring 
integral  equations. 


where 


P 

dP 

a + b P + C P^ 


o 


b a a/A 

o 

c a (y^  a/t  - v)/2  1. 

o o 


(1-14) 
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The  solution  of  Equation  (1-14)  depends  upon  the  sign  of 


If  q > 0,  one  obtains  an  inverse  tangent  solution  whereas  q < 0 yields 
an  inverse  hyperbolic  tangent.  However,  Equations  (1-12)  and  (1-13) 
imply  that  q < 0 in  order  for  the  restriction  to  be  satisfied.  There- 
fore, only  the  solution  for  q < 0 will  be  presented  here,  i.e.. 


A .1 

B tanh 

B 2-^  (t  - t^)  + tanh  (1/B) 
o 

- 1 

(1-15) 


Note  that,  because  the  slope  of  the  tanh  function  only  goes  to 
zero  as  t the  maximum  penetration  given  by  Equation  (1-12)  occurs 
theoretically  at  t However,  for  practical  purposes  P is  within 

1%  of  P^jax  argument  has  a value  of  only  2.65.  This  result  is 

analogous  to  that  obtained  by  hydrodynamic  codes  for  impact  problems, 
i.e.,  the  penetration  velocity  becomes  extremely  small  but  never 
reaches  zero. 


APPENDIX  II 


AXIAL  STEADY-STATE  SOLUTION  AND  ORDER-OF -MAGNITUDE  ANALYSIS 

In  this  case  we  assume  that  V.  s 0,  U « 0,  and  H > 0 and  Equation 
(6)  yields 


0 ■ [o  + (Vj  - U)^p^]  - - FjHn  2U 

(II-l) 

and  if  the  usual  jet  assumptions  of  negligible  strength  and  viscosity 
are  made,  one  obtains 


Vj  = (1  + Y)U  (II-2) 

where  y = v^Pj/p . This  is  the  well  known  shaped-charge  jet  expression^'**. 
Using  the  derinition  of  U * dP/dt  and  the  fact  that  V is  a constant 
(hence  not  given  by  Equation  8b)  one  obtains 


P = Vj(t  - t^)/(l  + Y).  (II-3) 

Equation  (1)  can  now  be  used  with  Equation  (II-3)  to  find  the  time  t^(2) 
when  the  penetrator  length  is  consumed,  i.e., 

t^(2)  - t^  = (1  + Y)  . (1 1-4) 

Substituting  Equation  (II-4)  into  (11-3)  >ields 


P 


max 


Vv 


in-53 


which  is  the  classical  jet  relationship.  Note  that  is  not 

defined  since  U is  assumed  to  be  a constant. 

Bernoulli's  equation  may  also  be  obtained,  without  recourse  to  the 
steady  state  assumption,  by  using  a nondimens ional  form  of  Equations 
(6) , (1)  and“(8b) , 


H r^^  U + H r^^  H/2  * r^^  [a  - a.^.] 


-U)2.^u2] 


(11-6) 
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and 


V “ ^ Cl  ■ 


1 ^ P - Ct  - t,) 

1 . (t  - t^)(l  - 


(11-7) 


where  the  *'bar*'  over  each  term  represents  a nondimer sional  value.  These 
equations  were  normalized  using  )Iq  as  the  characteristic  length,  Vq  as 
the  characteristic  velocity,  and  Iq/Vq  as  the  characteristic  time. 

d V p 

We  have  introduced  the  Reynolds  number,  Re_  ® -2 — 2 — 2.  ^ ^he  normalized 

— 2 — 2 
strengths  (the  Best  or  Metz  numbers)  o = o/p^  and  o.j,  = 

and  the  density  ratio  as  the  appropriate  nondimensional  groups. 

Now  we  define 

r.  = r./t  < < 1 or  r.  = 0*  (m)  where  m is  small,  H = 0(m), 
jjo  j 

Vj  = 0(1),  U = 0(1),  U is  of  0(1),  or  less,  and  for  H of 
0(m),  or  less,  we  have 


(Vj  ^U)2  - = 0 


(11-8) 


ignoring  terms  of  order  m or  less.  In  this  case  p.^,  /p  is  assumed  to 

Pt  p 

^ o “*2  2 

be  of  order  one,  or  in  general  U is  of  the  same  order  as  (V_  - U) 

% 1 

since  there  is  no  a priori  reason  to  assume  otherwise.  The  Best  numbers 
are  assumed  to  be  less  than  one,  i.e.,  ©(m)  and  the  Reynolds  number  is 

assimied  large,  i.e.,  of  9(-^). 

m 

The  equation  for  is  not  altered  by  the  order-of-magnitude 
analysis. 


*The  symbol  Q denotes  order~of -magnitude. 
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Now  for  continuous  jet, 


4 ^ 


% « ~ (V  - V 
o V ' o ro"^ 
o 


and 


Vt  T— 

U ■ -j-j — , Y ■ /pQ  via  the  order-of-inagnitude  analysis. 


(II-9) 


and 


> o 


Z + P 
o 


Vj  * — ^ — using  Equation  (II-9),  where  t^  * 

Since,  U = P,  the  above  equations  can  be  combined  to  yield 
/ dP  2o  " P 


dt  t(l  y)  * 


(II-IO) 


Integration  of  Equation  (II-IO),  for  P = 0 at  t * t^,  yields 


(II-ll) 


which  agrees  with  Allison-Vitali**  and  the  DSM^  continuous  jet  theory 
for  t < Tg,  (the  breakup  time). 


The  radial  flow  equations  (18)  - (23)  are  not  altered  by  the 
order~of -magnitude  analysis. 
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APPENDIX  III 


A.  Johnson *s  KE  Penetrator  Model  (Ref.  10) 

The  special  case  of  the  KE  penetrator  of  Reference  10,  which 
assumed  a constant  mass  projectile  with  a given  Initial  velocity,  Is 
considered  here.  Since  the  mass  Is  assumed  constant.  Equation  (1)  implies 
t = 0 or  U » Vj. 

In  this  case,  Equations  (6)  and  (9),  for  inviscid  flow,  yield 

+ H)  U = - Ot  - Pt  . (III-l) 

O 

Changing  the  independent  variable  from  t to  P,  since  U * P,  and 
neglecting  o^.  Equation  (III-l)  gives 


U - exp  (-  P/B  + H)) 


Cin-2) 


wnere  6 « Pq/p-j-  • Also, 
o 


and 


P - 6(«„ 


* H)  In 


o 

+ H)6 


(t  - 


^o5 


.l)  . 


(III-3) 


U 


1 + 


+ Hi 


(III-4) 


results.  Equation  (II 1-2)  illustrates  the  same  functional  form  as  given 
in  Reference  10.  It  is  interesting  to  note  that  Equation  (III-4)  is 
of  the  same  form  as  that  derived  in  Reference  4 for  the  penetration 
velocity  of  a particulated  shaped-charge  jet.  Furthermore,  both  the 
time  tjj(l)  (penetration  velocity  goes  to  zero  at  this  time)  and  the 
penetration  P go  towards  infinity  for  this  special  case. 
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Inclusion  of  the  strength  effects  in  these  equations  yields 


and 


P = + H)  tn 


(III-6) 


(t  ■ t ) 

where  9 = rj — rrr • In  this  case,  the  time  t (1)  is  given  by 

Po  + HJ  m 

the  following  equation. 


t (1)  - t 
■ m ' 


o 


(III-7) 


and  the  maximum  penetration  is  finite.  If  Equation  (III-7)  is  sub- 
stituted into  Equation  (III-6),  the  equation  for  maximum  penetration 
becomes 


Uo“  + 


0 . 

(III-8) 
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B.  Penetrator  with  Constant  Velocity 

It  is  also  of  interest  to  consider  a jet  without  a velocity  gradient; 
^ro  " Vq  a V,  or  a constant  velocity  rod.  In  this  case  Mjj  is  constant 
since  r j (t)  = r^  from  Equation  (3) , and  Equation  (6)  is  of  the  form 

U « A + BU  + CU^  (in-9) 


where  A 


. B 


2V 

H 


and  C = i 


2 

Letting  q = 4AC  - B , we  will  first  assiime  q > o,  and  in  this  case 
Equation  (III-9)  yields  the  following  solution, 


_2_ 


-1  2 C U + B -1  2 C + B 
tan tan  


(III-IO) 


When  q < 0,  which  is  usually  the  situation  since  the  jet  velocity  terms 
dominate  the  strength  and  viscosity  terms,  the  solution  to  Equation 
(III-9)  gives 


i-jriiuD  ,2CU+B 

tanh"^  ^ ^ - tanh"^  - ° 


(III-ll) 


Equation  (I I 1-9)  represents  an  analytical  derivation  of  the 
generalized  Poncelet  form  of  the  penetration  equations  which  has  been 
assumed  by  previous  investigators' 0 . 


It  is  noteworthy  that  when  the  jet  and  target  are  of  the  same 
density  (B  = 1)  Equation  (III-9)  yields 


(111-12) 


Thus,  3=1  may  be  easily  treated  as  a special  case. 
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For  e * 1 either  Equation  (III-IO)  or  Equation  (III-ll)  can  be 

illustrate  the  solution  for  the  case  of  Equation 
(I**-!!)  which  yields  the  following  penetration  velocity, 


2 C U » B 


(X  - tanh  0)/l  - X tanh  0) 


(III-13) 


and  the  corresponding  penetration  is, 


where 


2 C P + B (t  - t^)  = - 2 in  {cosh  6 - X sinh  0} 


0 = 


^ (t  - 


2 C U + B 


and  X = 


-q 


(III-14) 


The  length  of  the  penetrator  becomes. 


^ - V(t  - t^)  - 


B(t 


"o) 


2C 


- ~ tn  {cosh  0 - X sinh  0}.  (III-15) 


Inspection  of  Equation  (III-13)  shows  that  U remains  finite  as 
t -»■  <*>  and  hence  fne  penetration  given  by  Equation  (II I -14)  -*■  ».  jn  this 
case,  the  time  t^(i)  ® based  upon  the  penetration  velocity  U,  whereas 

^m(2j  IS  finite  since  2,(t)  -»■  o at  some  finite  time. 

9 and  10  compare  the  analytical  results  with  the  experimental 
data  of  Tate  for  aliiminum  and  duraluminum  rods  impacting  polyethylene 
targets.  Figure  9 presents  the  penetrator  erosion  as  a function  of 
time  where  the  analytical  model  includes  either  density  effects  (v) • 
density  ^d  strength  effects  (y,  a);  or  density,  strength  and  viscosity 
(Y.  o,  V)  effects  . The  analytical  results  bracket  the  experimental 
data  and  the  agreement  is  good.  In  Figure  10,  penetration  versus  time 
results  are  presented.  The  analytical  models  agree  with  the  experimental 
data  and  the  predicted  final -penetrations  for  density  and  strength  effects 
agree  very  well  for  both  the  hard  and  soft  aluminum-penetrators.i"*'  Note 
that  the  actual  final-penetration  for  both  penetrators  exceeds  that 
predicted  by  the  hydrodynamic  theory. 


■•’2^6  strength  of  alvminm}^  is  60  MPa  and  the  strength  of  duraluminum^^ 
^s  300  MPa.  Polyethylene  was  assumed  to  have  a density  of  0.92  gm/am^ 
and  zero  strength  . The  H-nematio  viscosity  of  polyethylene  was 
assumed  to  be  (0.5  mVeecJ^  i.e.,  close  to  that  of  lead,  since  no 
polyethylene  viscosity  data  could  he  found. 

■f'+rhe  shock  Hugoniot  velocity  of  1.1  km/s  was  used  as  the  initial 
penetration~velocity  U^. 


48 


APPENDIX  IV 


A CONSTANT  LENGTH  (NON-ERODING)  ROD  PENETRATOR 


This  Appendix  considers  another  special  case,  namely  that  the  rod 
length  is  constant  or  from  Equation  (1)  £.  » o or  ii.  = Aq  and  the  penetrator 
velocity  is  equal  to  the  penetration  velocity  or  U * VcG(i)*  Hence  the 
rod  mass  is  constant  (since  rj  = r^). 

Equations  (6)  and  (9)  may  be  combined  to  yield 

U (t^  + K)  = - [a  + bU  + CU^],  (IV-1) 

where  the  viscous  effects  for  region  I are  assumed  to  be  much  less  than 


°T  2H  u 

the  viscous  effects  for  region  II  and  where  a = — , b = j » 

*^o  p r 

Pt  ‘'o  o 

o 2 * 

c = - — and  q * 4ac  - b . Equation  (IV-1)  has  the  same  functional 
^0 

form  as  the  generalized  Poncelet  equation  which  has  been  used  exten- 
sively in  penetration  mechanics  studies  to  curve-fit  the  experimental 
dal  -' . 


Equations  (6)  and  (9)  are  combined  since  regions  I and  II  are  no 
lu'i7.ir  distinct  when  U ■ Vj.  These  two  equations  provide  one  equation 
in  'u,. 


Equation  (IV-1)  can  be  solved  for  q > o or  q < o as  in  Appendix  III. 
For  q < o,  i.e.,  the  viscous  terms  dominate,  the  penetration  velocity 
is 


2 Me  + b 


A + tanh  6 
1 + A tanh  6 


(IV-2) 


(t  - 1 ) 
® = 2 . H) 

t (1)  for  maximum  pent 


2 U c + b 
and  A = 2 

iion  to  be 


When  U = o we  determine 


■*For  inoiaoid  flow^  the  reautta  of  Appendix  IV  are  identical  to  thoae 
of  Appendix  III-A,  Appendix  III-A  was  preaented  aepaxvLtely  in  order 
to  Qompare  the  Appendix  III-A  reautta  with  Ref,  (10), 


49 


and  the  penetration  is  given  oy 


2P  c + b (t  - t^)  » 2 + H)  An  {cosh  0 + A sinh  0).  (IV-4) 

The  maximum  penetration  is  obtained  by  substituting  Equation  (IV- 3)  into 
Equation  (IV-4) . 

When  q > o,  i.e.,  the  strength  and  density  terms  dominate,  we  have 


2 c U + b - y^A  - 

1 + A tan  6 


(IV-5) 


where  now  0 = 


(t  - t^) 

TIFT-RT 


2 U c + b 


and  A » 


. The  maximum  time  to 


penetrate  tj|j(l)  becomes, 


t (1)  - t = 
m^  o 


2 * w -I  r “o  1 


CIV-6) 


This  equation  shows  that  the  time  to  reach  maximum  penetration  increases 
linearly  with  penetrator  length,  and  decreases  in  a non-linear  fashion 
•with  specific  target  strength  (oj/pj  ) . Equation  (IV-6)  also  shows  that 
t^(l)  approaches  a limiting  value  as^U^  becomes  large.  The  penetration 
is 


, (Jl  + H) 

P = -^(t-t)  + —2 Hi,  {cos  e + a sin  6}.  (IV-7) 

M ^ 0 C 


Inserting  Equation  (IV-6)  into  Equation  (IV-7)  yields  the  general 
expression  for  maximum  penetration. 


The  analytical  results  of  this  Appendix  are  compared  with  the 
results  of  Reference  14  for  a steel  bullet  impacting  dry  sand.  The 
penetrator  may  be  realistically  assumed  to  be  non-eroding.  In  this 
case.  Figure  11  shows  good  agreement  between  the  predicted  and  experi- 
mental penetration  versus  time,  penetration  velocity  versus  time,  maxi- 
mum penetration  depth,  and  time  to  reach  final  penetration  depth. 


"^The  Band  density  is  1.68  gm/am^,  the  etrenqth  of  sand  ia  10^  Pa^  and 
U = 0.6?i  km/aeo  waa  taken  for  the  initial  Hugoniot  velocity.  The 
aenaity  of  ateel  ia  7.8  gm/cnr  and  inviaoid  flow  waa  aaavmd. 


50 


m 


Figure  12  presents  the  penetratlon-tlne  history  for  a non-eroding 
copper  penetrator  impacting  aluminiun,  zinc,  steel  and  lead  targets. 

Ail(l)  or  Fe(l)  denotes  a lower  limit  on  the  target  strength  whereas 
At (2)  or  Fe(2)  denotes  the  higher  strength  term  for  the  specified  target 
material.  Standard  handbook  values  for  the  strength  terms  were  utilized. 
Note  that  the  trend  of  the  data  agrees  with  the  experimental  results 
from  Figure  7 (after  the  normalization  is  removed) . However,  the  non- 
eroding penetrator  model  overpredicts  the  actual  penetration-time 
performance. 

Figure  13  compares  the  analytical  results  with  the  data  of  Reference 
12  for  a copper  rod  impacting  an  annealed  copper  target.  Note  that  the 
analytical  results  agree  closely  with  the  various  experimental  measure- 
ments.'^ 

Finally,  Table  IV-1  summarizes  the  maximum  predicted  penetration 
values  for  a copper  penetrator  impacting  the  noted  targets.  The 
analytical  results  follow  from  Appendices  II,  III  and  IV  and  the 
experimental  data  is  from  Figure  7.  The  subscript  L denotes  the  lower 
limit  of  the  variable  in  question  and  a subscript  U denotes  the  upper 
limit.  The  standoff  is  denoted  by  2^  and  represents  the  maximum 

length  (TgCVg  - Vj,q))  of  the  fully  stretched  jet.  In  spite  of  the 
various  assumptions  associated  with  the  different  models,  the  analyti- 
cal data  brackets  the  experimental  results  and  is  of  the  same  order  of 
magnitude.  Thus,  the  authors  are  optimistic  that  numerical  solutions, 
which  will  remove  the  limiting  assumptions,  will  accurately  predict 
shaped-charge  jet  penetration  and  the  resulting  target  response. 


'^Irw'iaoid  floD  waa  aaawred  and  the  initial  Hugoniot  velocity  was  used 
(U^  - 0.57  hn/aeo). 
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r. 

3 


CG 


r = r(t  ) 


height  of  Interaction  region 
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moving  mass  of  eroded  penetrator  per  unit  radial  segment 

penetration 

maximum  penetration 
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moving  target  material 

radial  location  of  the  center-of-mass  target  material 

Reynolds  number  = p V S,  /y 

0 0 

radius  of  the  hole  profile  at  some  depth  P 
penetrator  radius 

radius  to  center-of-mass  of  penetrator  material  within 
the  interaction  region 

initial  penetrator  radius 

inside  radius  of  the  eroded  penetrator  material 

time  for  penetrator  to  reach  target 

time 
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